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Abstract 

Rank metric codes and constant-dimension codes (CDCs) have been considered for error control in 
random network coding. Since decoder errors are more detrimental to system performance than decoder 
failures, in this paper we investigate the decoder error probability (DEP) of bounded distance decoders 
(BDDs) for rank metric codes and CDCs. For rank metric codes, we consider a channel motivated by 
network coding, where errors with the same row space are equiprobable. Over such channels, we establish 
upper bounds on the DEPs of BDDs, determine the exact DEP of BDDs for maximum rank distance 
(MRD) codes, and show that MRD codes have the greatest DEPs up to a scalar. To evaluate the DEPs of 
BDDs for CDCs, we first estabUsh some fundamental geometric properties of the projective space. Using 
these geometric properties, we then consider BDDs in both subspace and injection metrics and derive 
analytical expressions of their DEPs for CDCs, over a synmietric operator channel, as functions of their 
distance distributions. Finally, we focus on CDCs obtained by hfting rank metric codes and establish 
two important results: First, we derive asymptotically tight upper bounds on the DEPs of BDDs in both 
metrics; Second, we show that the DEPs for KK codes are the greatest up to a scalar among all CDCs 
obtained by lifting rank metric codes. 

I. Introduction 

Two important classes of codes for error correction in random network coding are constant-dimension 
codes (CDCs) and rank metric codes. The first class of codes are CDCs, which refer to codes defined in 
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Grassmannians associated with the vector space over a finite field. They correct the errors and erasures 
which typically occur on the network by using the subspace metric ||TJ, while using the injection metric 
they correct errors on adversarial channels ||2|. Many bounds and constructions were proposed in the 
literature (see, for example, ||T|, Q, Q), and CDCs were shown to be asymptotically optimal subspace 
codes |5|. The second class of codes are rank metric codes ||6|-|[8|, which are related to CDCs in 
more than one way. First, nearly optimal CDCs can be constructed by lifting rank metric codes |9|, a 
transformation that preserves the distance. Second, it was recently shown that CDCs are closely related 



to constant-rank codes |10|. More importantly, error correction in coherent and noncoherent random 
network coding can be solved from a rank metric perspective |[2|, |[9|. There has been a steady stream 
of work on rank metric codes, in particular the maximum cardinality of a code with a given minimum 
rank distance was determined in |[6|-|[8|. We refer to codes with maximum cardinality as maximum rank 
distance (MRD) codes, and the class of linear MRD codes proposed independently in |[6|-|[8| as Gabidulin 
codes henceforth. 

One critical aspect that has received little attention in the literature is the error performance of 
bounded distance decoders (BDDs) for rank metric codes and CDCs. Given a received word, a BDD 
either declares a failure or finds a codeword within a predetermined radius of the received word. In the 
latter case, when the codeword produced by the BDD is not the sent codeword, a decoder error occurs. In 
many applications and especially in network coding, a decoder error is more detrimental than a decoding 
failure, and the decoder error probability (DEP) is a crucial parameter of BDDs for the codes used in 
the transmission. 

In this paper, we investigate the DEP of BDDs for rank metric codes and CDCs on channels motivated 
by error control in random network coding. The main contributions of this paper follow. 

• We first consider a channel on matrices where all additive errors with the same row (or column) 
space are equiprobable. For any rank metric code over the equal row (or column) space channel, we 
derive upper bounds on the DEP of bounded rank distance decoders. Our results show that the DEP 
of BDDs for any rank metric code decreases exponentially with the square of its error correction 
capability. For MRD codes over the equal row (or column) space channel, we derive the exact DEP 
of BDDs and show that MRD codes have the highest DEP up to a scalar. 

• We then consider operator channels on subspaces, more precisely symmetric operator channels, 
where all outputs with the same dimension and at the same subspace (or injection) distance from 
the input are equiprobable. Using geometric properties of balls with subspace radii, we derive the 
DEPs of BDDs for arbitrary CDCs over a symmetric operator channel. The analytical expressions 
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for both metrics ultimately depend on the distance distributions of CDCs. 

• Finally, we focus on CDCs obtained by lifting rank metric codes since their distance distributions 
are known, and obtain two important results. First, we obtain asymptotically tight upper bounds on 
the DEPs of BDDs in both metrics; the upper bounds decrease exponentially with the square of the 
decoding radius. Second, we show that the DEPs for KK codes, which are nearly optimal CDCs yj 
and can be obtained by lifting Gabidulin codes |[9|, are the highest up to a scalar among all CDCs 
obtained by lifting rank metric codes. 

Our work on the error performance for CDCs is novel to the best of our knowledge. Our work on 



the error performance for rank metric codes differs from the work in 1 1 1 1 in several aspects, and is a 
nontrivial extension of our previous work [12]. The error performance analysis in |11| was aimed at 
two-dimensional errors in data storage equipments and as such, assumes different channel models and 
considers decoder errors and decoder failures together. Our results in this paper are more general in terms 
of both the channel model and underlying codes than those in 1 12], and they required a more fundamental 
geometric approach and the use of novel techniques such as constant-rank codes pO} . This more general 
investigation of decoder error performance is important, since explicit construction of optimal constant- 
dimension codes for arbitrary parameters is unknown and suboptimal codes are sometimes favored due 
to complexity or error performance. 

Our work on the error performance for rank metric codes parallels some previous works for Hamming 



metric codes. In |13|, an upper bound on the DEP of a bounded Hamming distance decoder for linear 
MDS codes over channels where all errors with the same Hamming weight are equiprobable was derived. 
This work was refined in y_4J, where the exact DEP for linear MDS codes was determined under the 
same setting. In | |T5[ , the results in |fT3l were extended to more general channels and to any linear code. 
More precisely, p3| introduces error-value symmetric channels, where all errors with the same support 



are equiprobable, thus taking bursty channels into account. Our results on the DEP for rank metric codes 



over the equal row (or column) space channel parallel those in 1 15 1. We would like to emphasize that the 
work in this paper strongly differs from the previous work reviewed above. First, while readily available 
combinatorial results can be used to derive the results for Hamming metric codes (for example, the results 



in 1 14 1 are based on inclusion-exclusion principle), their counterparts for rank metric codes have to be 
established explicitly. Also, our results for rank metric codes and CDCs are based on completely different 
approaches from those used for Hamming metric codes. 

The rest of the paper is organized as follows. Section [IT] reviews some necessary background. In 



Section pll we investigate the DEPs for rank metric codes. In Section IV we investigate the DEPs for 
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CDCs in general and liftings of rank metric codes in particular. 



II. Preliminaries 



A. Rank Metric Codes 



def 



The rank distance between two matrices in GF(g)™^" is defined as ^^(X, Y) = rk(X — Y). The 
number of matrices at distance u from a given matrix is denoted as A'r(u) = where 
a{m,u) = for m < 0, a(m,0) = 1, and a(m, u) = nr=o^(^'" ~ 1^) for m > and u > 1, and 



a{n,u) 
a{u,u) 



is referred to as the Gaussian binomial. For all < r < n, we have 



r(n-r) ^ 



< 



^-1 r(n-r) 



(1) 



where Kg = YlJ^iil — q~^) 1 12 1. The volume of a ball with rank radius t in GF(g)™^" is denoted as 
V^{t) = X]s=0"^r(^)- "^^'^ intersection number J^{u,s,d), defined as the volume of the intersection of 
two spheres with radii u and s and distance d between their centers, was derived in p6|. In particular, 



the intersection numbers satisfy jR(t, d — t,d) = q^^'^ [^] for all < t < d < min{n, m}, and 1 17 1: 

N^{d)Mu,s,d) = N^{u)A{d,s,u) (2) 

n 

Y,Mu,s,d) = N^{s). (3) 



u=0 



A rank metric code can be viewed as a subset of GF((7)™^", where the minimum rank distance of a 
code is simply the minimum distance over all pairs of distinct codewords. The maximum cardinality of 
a rank metric code in GF((7)™^" with minimum rank distance d is min{g'"("^'^+^\ g"^'""'^"'"^)} |[6|-|[8|. 
We refer to codes with maximum cardinality as maximum rank distance (MRD) codes. The number of 
codewords at rank distance r from a given codeword in an MRD code in GF(g)'"^" (n < m) with 
minimum rank distance d was determined in |[6|, ||7| and is denoted as M{d,r). In particular, we have 
M(d,d)= [a(g--l). 



A constant-rank code is a rank metric code whose codewords have the same rank pO] |. The maximum 
cardinality of a constant-rank code in GF{q)^^^ with minimum rank distance d and rank r, denoted as 



Af^{q, m, n, d, r), satisfies A^{q, m, n, d, r) = Af^{q, n, m, d, r) and, forn < m and d < r 1 10 Proposition 
8], 



A^{q,m,n,d,r) < 



a{m, r — d + 1). 



(4) 



'For completeness, it should be denoted as NR{q,m,n,u). When there is no ambiguity about the underlying field and the 
size of matrices, we use Ns.{u) for simplicity. To simplify the notation in this paper, similar simplifications have been made to 
other symbols. 



April 1, 2010 



DRAFT 



5 



B. Constant-Dimension Codes 

We refer to the set of all subspaces of GF(g)" with dimension r as the Grassmannian and denote it 
as Er{q,n). We also refer to the projective space [J^=QEr{q,n) as E{q,n). For U,V £ E{q,n), their 
intersection U CiV is also a subspace in E{q,n), and we denote the smallest subspace containing the 
union of U and V as U + V. The subspace metric (I] and the injection metric |[2|, respectively defined 
as 

ds{U,V) = dim(C/ + y) - dim(C/ n y) 

d,{U, V) = ^ds{U, V) + ^\ dim(C/) - dim(y)| = max{dim(C/), dim(y)} - dim(C/ n V) 

are both metrics over E{q,n). 

The number of subspaces with dimension s at subspace distance d from a subspace with dimension r 
(0 < r < n), denoted as Ns{r,s,d), is [^] [^1^] when u = ^-+^^ is an integer, or otherwise 

Q. The number of subspaces with dimension s at injection distance d from a subspace with dimension 
r is given by Ni{r, s, d) = Ns{r, s, 2d — \r — s\) [5]. In the Grassmannian Er{q, n), the intersection of 
any two spheres of injection radius u and s with injection distance d between their centers, referred to 
as the intersection number Jc(r, n, s, d), was determined in 

A subset of Er{q,n) is called a constant-dimension code (CDC). Since dfi{U,V) = dfi2di{U,V) if 
U and V have the same dimension, the minimum subspace distance of a CDC is equal to twice its 
minimum injection distance. CDCs are related to rank metric codes through constant-rank codes or 
through the lifting operation jol, described below. The lifting of C G GF(g)''^('^-'') is defined as /(C) = 
R{Ir\C) G Er{q, n), where I,, is the r x r identity matrix and R denotes the row space of a matrix. For 
all C,D G GF(g)'^^("-'^), we have (is(/(C), /(D)) = 2d,(/(C), /(D)) = 2^(0, D) f^. Therefore, the 
injection distance distribution of the lifting is equal to the rank distance distribution of the original code. 
Liftings of MRD codes were introduced in ||9|, and we refer to them as KK codes. 

III. Decoder error probability of rank metric codes 

In this section, we investigate the DEP of bounded rank distance decoders for rank metric codes. All 
distances are rank distances in this section. 

We assume the following scenario, where an adversary injects linearly independent packets mali- 
ciously on the network, using some knowledge about the transmitted packets or the protocol used. Hence 
the adversary may choose and inject linearly independent packets so as to corrupt the transmitted message 
more effectively than others (see the example in [2]). These packets undergo linear- combinations through 
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the network, and result into an additive error at the receiver whose rank depends on the number of packets 
injected. Due to the vector-space preserving property of linear network coding, the row space spanned by 
these packets remains unchanged through the Unear combinations operated at intermediate nodes. Hence, 
the additive error at the receiver can take any value provided its row space is fixed. This leads to our 
model of additive errors, where all errors with the same row space are equiprobable. Also, because the 
rank is preserved by transposition, we also consider channels where errors with the same column space 
are equiprobable. A channel on GF(g)™^" is hence said to be equal row (column) space if errors are 
additive and the errors with the same row (column) space are equiprobable. 

Although equal row and column space channels were proposed to model error correction in random 
linear network coding as described in Section I, it is remarkable that these channels may also be used to 
model other appUcations of rank metric codes. Rank metric codes can be used for the correction of two- 
dimensional errors |[8|, 1 18 1 (i.e., errors confined to a certain number of rows and columns) in storage 
equipments. Hence, our model encompasses the case of two-dimensional errors, where some rows or 
columns are more likely to be in error than others. 

Note that using a code C over an equal row space channel on GF(g)'"^" is equivalent to using its 
transpose code over a column space channel on GF(g)"^™. Since the transpose operation preserves the 
distance, we only study equal row space channels. In this case, the DEP depends on the sent codeword 
C, the row space U G Eu{q,n) of the error, and the error correction capability t of the code. We hence 
derive a bound on Pr(C, f7, t) for all codes, and we also obtain the exact value of the DEP for linear 
MRD codes when n < m. Since the distance distribution of an MRD code does not depend on the 
codeword, let us denote the DEP of a BDD for an MRD code in GF(g)™^'^ with minimum distance d 

as P^MwiU, t). 

Let P^{C,t) the DEP of a bounded distance decoder with radius t for the code C. We have 
P^{C,t) = Y, PR(C,C/,t)P{C,[/} <maxPR(C,C/,t), 

CeC U&E{q,n) 

where P{C,U} is the probability that the transmitted codeword is C and that the channel error has row 
space U. Therefore, although the probability Pr{C, U, t) is conditional, it can be used as an upper bound 
on the unconditional DEP. The upper bounds on Pr{C, U, t) that we shall derive can hence be applied 
for the overall DEP of any rank metric code. 

For any R G GF{q)"^^^, we denote the number of matrices with row space W and at rank distance 
s from R as g^CW, 8,11). We prove below that gn{W, 8,11) depends on R only through its row space. 

Lemma 1: For all R, S G GF(g)™^" with the same row space U, g^^iW, s, R) = 5'r(VF, s, S). 



April 1, 2010 



DRAFT 



7 



Proof: Suppose X € GF((7)™^" has row space W and satisfies rk(X — R) = s. Expressing 
S = AR where A G GF(q)™^'" has full rank, the matrix Y = AX has row space W and satisfies 
rk(Y - S) = rk(X - R) = s. Thus g^{W,s,S) > g^{W,s,K). Using R = A^^S, we show that 
9k{W, s, S) < g^{W, s, R); hence g^{W, s, R) = g^{W, s, S). ■ 
Since 5r(VF, s, R) depends on R only through its row space U, we denote it as 5'r(VF, s, U) henceforth. 

def 

The DEP is derived in terms of the distribution of codewords according to their row space Ay^iC) = 
|{D G C : i?(D -C) = W}\ for all W £ E{q,n) and the g^iW, s, U) constants in Proposition [T] below. 

Proposition 1: Assuming a codeword C € C, a code in GF(g)'"^" with minimum distance d, is sent 
over an equal row space channel and the channel error has row space U € Eu{q, n), the DEP of a BDD 
with decoding radius t = [^^J is when u < d — t, and when u > d — t satisfies 

^ n t 

P^iC,U,t) = E Aw{C)J29r{U,s,W) (5) 

a[m,u) ^-^ ^-^ ^-^ 



1 " r 1 * 

< , , ^ a(m,u;-d+l)^ Jr(u, s,?/;) (6) 



w 



^ 2 -t(m-n+t) ^j^gj^ d = 2t + 1, 

< < ^ (V) 

^-2^-t(m-n+t)-m ^j^^j^ d = 2t + 2. 



Furthermore, if n < m and C is a linear MRD code, then the DEP is given by 

t 



PkmkAU, t) = — ^ ^ M{d, w) J2 Mu, s, w). (8) 



w=d s=0 



The proof of Proposition [T] is given in Appendix [A] We remark that the upper bound in ([6]) actually 
depends on u only, not U, while the bound in (|7]) is very general: it does not depend on the transmitted 
codeword C and the rank of error u, and applies to any rank metric code. In fact, applying ([T]) to linear 



MRD codes leads to |12 Proposition 6]. 

We now show that in all nontrivial cases MRD codes have the greatest DEP among all rank metric 
codes up to a scalar. 

Corollary 1: Let C be any rank metric code in GF(q)™^" (n < m) with minimum distance d and 
let C G C. Then if g > 2, n < m, or ci 7^ m - 1, Pr(C, U,t) < HqPf^uRDiU,t), where H2 = 3.5 and 
Hq='^forq> 2. 



Proof: By 1 10 Proposition 15], we have HqM{d, r) > Aji{d, r) for n > r > d, provided that q > 2, 

n < m, ox d ^ m - I. Hence HqP^Mm{U, t) > Y2,=d ^^i^, w) Y!s=q M'^^ w) > Pr(C, U, t). 
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We remark that Corollary [T]does not hold only if q = 2, n = m, and d = m — l = n — 1, which is 
a trivial case. Corollary [T] indicates that there is a tradeoff between decoder capability (or radius) and 
decoder error probability. Given m, n, and the dimension, MRD codes achieve the greatest minimum 
rank distance and hence have the greatest decoder capability. However, they also have the greatest DEP 
among all rank metric codes up to a scalar. 

We would like to emphasize our choice to consider the DEP of any rank metric code, and not only 
MRD codes for which an efficient BDD is known. In our scenario, the DEP corresponds to a security 
criterion: it is the probability that the adversary manages to thwart the communication. A low DEP is 
hence equivalent to a high protection against the attacks by the adversary. Since MRD codes have the 
highest DEP up to a scalar, they may not provide a strong enough level of security. It is hence reasonable 
in appUcations to consider another class of suboptimal codes in order to increase the level of security. 
Although decoding algorithms have been proposed for MRD codes only, we show below that suboptimal 
codes could also be used with a low-complexity decoder. Such codes can be easily designed by taking 
Gabidulin codes over a subspace. Also, in his seminal paper, Gabidulin introduced a whole class of codes, 
referred to as g-cyclic codes, which are analogues of cyclic codes. Amongst these codes, one can find the 
analogues of BCH codes and in particular, the analogues of Reed-Solomon codes which are Gabidulin 
codes. Since many algorithms analogous to the ones for Reed-Solomon codes have been proposed for 
Gabidulin codes, the analogues of BCH codes are likely to also have low-complexity decoders. Using 
such suboptimal codes may yield a much lower DEP, and hence a higher level of security. 

In the particular case of a channel where all errors with the same rank are equiprobable (referred to as 
rank symmetric in 1 12]), the DEP only depends on the rank u of the error. It can be easily shown from Q 
that it is related to the distance distribution y4«,(C) = |{D £ C : rk(D-C) = w}\ = J2weE„(q n) ^w{C) 
by 

^ n t 

P^{C,u,t) = -—y2AMy2Mu,s,w). (9) 

w=d s=0 

The result in Q is interesting, as the distance distribution of a code has been widely studied, and the 
distance distribution is usually known for a given code, while the distribution of codewords according to 
their row space appearing in (|5]l has received much less attention. 

In many applications, the probability that the received matrix is at distance u from the sent codeword 
decreases rapidly with u, and hence the overall DEP can be approximated by Pj^{C,d — t,t). We note 
that d — t is the smallest value for u that will lead to a decoder error. Furthermore, when u = d — t, the 
channel output can only be decoded to codewords that are at distance d from the sent codeword. 
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Proposition 2: Assume a codeword C G C, a code in GF(g)™^" with minimum distance d, is sent 
over a rank symmetric channel and the channel output is at distance u = d — t from C, the DEP of a 
BDD with decoding radius t = [^^J is given by 

P,{C,d-t,t)=q'^'-' \ p A,iC). (10) 

In particular, the DEP for an MRD code satisfies 

f K when d = 2t + l, 

PKM^{d-t,t)> { (11) 
[ when d = 2t + 2. 

Proof: Since the channel output can be decoded to only codewords at distance d from the sent 

codeword C, (6l reduces to Pt,{C,d-t,t) = ^^^^^^^^^(C), which gives (jioj). For an MRD code, (8 1 

and the bounds on the Gaussian binomial in ([T]) yield ( [TT| . ■ 

When u = d — t, Proposition |2] above not only provides the DEP for any code, but also shows that 

the upper bound on the DEP for MRD codes in ^ is tight up to a scalar since the upper bound in ^ 



and the approximation in (111 differ by only a scalar which tends to 1 wen q increases. 



IV. DEP FOR CDCS OVER A SYMMETRIC OPERATOR CHANNEL 

A. Further Properties of Balls with Subspace Radii 

Properties of balls in E{q, n) with subspace or injection radii were investigated in Q. In this section, 
we determine further properties of such balls, which will be instrumental in our analysis of DEP for 
CDCs. 

We study the properties of balls with subspace radii only, as properties of balls with injection radii 
will not be useful to our derivation of the DEP. Proposition |3] below shows that the intersection of two 
spheres in the projective space only depends on the radii of the spheres, the distance between between 
the centers, and the dimensions of the centers. 

Proposition 3: For all A G Ea{q,n),B € Ei,{q,n) with d^{A, B) = w, the number of subspaces 
C G Ec{q, n) such that d^{A, C) = u and (is(-B, C) = s only depends on u, s, w, a, b, and c. It is hence 
denoted as Js(n, s, w; a, b, c). 

Proof: Let t = dim(A H B) = ^{a + b — w) and vq, vi, . . . j^Ta+b-t-i be linearly independent 
vectors such that vq, vi, . . . , vt_i e AnB, vj, vt+i, . . . , v^-i G A, and v^, Vf+i, . . .,Va+b-t-i e B. 
The matrix V G GF(g)('^+^-*)^" whose rows are given by the VjS has full rank, therefore there exists a 
nonsingular matrix X G GF(g)'^^" such that VX = h+a-t- Then AX. = {uX : u G A} is the span of 
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the unit vectors eg, ei, . . . , e^-i, denoted as /„ and 5X is the span of eo, ei, . . . , e(_i, e^, . . . , Ba+b-t-i, 
denoted as Jb- Therefore, ds{Ia, Jb) = ds{A, B) = w and 

\{C E E,{q,n)\d,{A,C) = u,d,{B,C) = s}\ 

= |{CX E Ec{q,n)\dsiAX,CX) = u,ds{BX,CX) = s}\ 

= \{Y E E,{q,n)\ds{Ia,Y) = u,dsiJb,Y) = 

which is a function of u, s, w, a, b, and c. ■ 
The intersection of spheres with subspace radii is related to their volume in Corollary |2] below. 
Corollary 2: For all parameter values, 

Ns{a,b,'w)Js{u,s,w;a,b,c) = Ns{a,c,u)Js{w,s,u;a,c,b), (12) 

n 

'^Js{u,s,w;a,b,c) = Ns{b,c,s). (13) 

u=0 

Proof: Let A E Ea{q, n). By counting the number of pairs of subspaces {B, C) such that dim(i?) = 
b, dim(C) = c, ds(j4,S) = w, ds,{A,C) = u, and ds{B,C) = s in two different ways, we obtain ( [T2) ). 
Next, let B E Ef,{q,n) be fixed; we denote ds{A,B) by w. As Y^2,=o "^si^, s,w; a,b, c) is the number 
of subspaces C E Ec{q,'n) such that d^{B,C) = s, this sum equals Ng{b,c,s). ■ 
Although the value of Js{u, s,w;a,b,c) is unknown in general, we remark that for a = 6 = c, 
Js{2u, 2s, 2w; a, a, a) = Jc(a, u, s, w), the intersection number for Ea{q, n). We also determine its value 
when w = u + s below. 

Proposition 4: Js{u, s,u + s; a, b,c) = when u > min{a + c, a + 26 — c}, s > min{6 + c, 2a + b — c}, 
or u + s > min{a + b, n}. Otherwise, we have 



Js(n, s,u + s; a, b, c) 
The proof of Proposition [4] is given in Appendix |B 



r a—b+u+s^ 




rb—a+u+si 


2 




2 


c—b+s 




c—a+u 


L 2 J 




L 2 -1 



(14) 



B. DEP of Bounded Subspace Distance Decoders for CDCs 

We study the DEP of a bounded distance decoder for a CDC in Er{q,n) over a symmetric operator 
channel, defined below. An operator channel is a channel where the inputs and outputs are subspaces in 
E{q,n). As assumed in the channel may erase some dimensions of the transmitted subspace as well 
as inject some erroneous dimensions. We refer to these as erasures and errors, respectively. We say an 
operator channel is symmetric when all outputs corresponding to e errors and p erasures are equiprobable. 
If the input has dimension r and e enors and p erasures occur, the output has dimension v = r + e — p 



April 1, 2010 



DRAFT 



11 



and is at subspace distance u = e + p and at injection distance = max{e, p} from the input. Therefore, 
an operator channel is symmetric if and only if all outputs with the same dimension and at the same 
subspace (or injection) distance from the input are equiprobable. 

Since the distance properties of a code C in Er{q,n) are equal to those of the code in En-riq^n) 



We first study the DEP of a bounded subspace distance decoder. We remark that u + v — r = 2e 
should be an even integer, otherwise there are no subspaces with dimension v and at distance u from the 
transmitted subspace, that is, N^{r,v,u) = 0. 

For a CDC C C Er{q,n) with minimum subspace distance 2d, suppose a codeword C is transmitted 
over a symmetric operator channel. The output of a BDD with decoding radius d— 1 and its DEP depend on 
both e and p, or equivalently on both u and v. First, if u < d—1, that is, the channel output is at subspace 
distance up to d — 1 from the sent codeword C, the BDD will produce the sent codeword C. When u = d, 
the channel output is beyond the decoding radius of any codeword, and hence the BDD will produce a 
failure. When u> d+1 and \ v — r\ > d—1, since the distance between the received subspace and the code 
is no less than r|, the channel output is beyond the decoding radius of any codeword and the BDD will 
produce a failure. In all the cases above, the DEP is zero. When u > d + 1 and r — d + 1 < v < r + d—1, 
a decoder error is possible, and we determine the DEP based on the distance distribution of the code 
with respect to C, denoted as A^{C) = \{D £ C : di{D, C) = w}\ = \{D £ C : ds{D, C) = 2w}\. 

Proposition 5: Assuming a codeword C of a CDC in Er{q, n) with minimum subspace distance 2d is 
sent over a symmetric operator channel and that the received subspace has dimension v and is at subspace 
distance u from C, the DEP of a bounded subspace distance decoder with decoding radius d — 1 is given 



when A''s(r, v,u) > {), u> d + 1, and \v — r\ < d — 1; and Ps{C, u,v,d—l) = otherwise. Furthermore, 
if the CDC is the lifting of a rank metric code, then 



consisting of the dual subspace of each codeword in C, we assume 2r < n as in yj. 



by 




r 



d-1 



(15) 



Ps(/(C),n,i-,d-l) 




(16) 




< < 



(17) 



where Lq = 
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The distinction ond— 1 + r — vin the upper bound in ( [TT] ) is explained as follows. For d — I + r — v 
even, the largest subspace distance of a decodable subspace is exactly d — 1; however, it is only d — 2 
when d— 1 + r — vis odd. Therefore, for v such that d — l + r — vis odd, the code uses a bounded 
distance decoder with true decoding radius d — 2 instead of d — 1, leading to a smaller DEP. 

Note that the bound in ( [TT] ) is very general, as it does not depend on the transmitted codeword or 
the distance of the received subspace to the sent codeword. We also remark that the exponent in the 



bound in (17) becomes zero when either v = r — d + lorr='^ and v = | + d — 1. The case 
where v = r — d + 1 can be explained as follows. The decoding region around each codeword consists 
of decodable subspaces with dimension v = r — d + 1, and hence is the intersection of the sphere of 
radius d — 1 around a codeword and the Grassmannian £'r_d+i(g', n). The volume of each decoding 
region is N^{r,r — d + 1, d — 1) = [/l^], and for the lifting of an MRD code, the disjoint union of 
the decoding regions has ^^'^^j > q{r-d+i){n-r+d-i) subspaces by (jT|. On the other hand, 

\Er-d+iiq,n)\ < q(r-d+i)(n-r+d-i) hcncc the decoding regions form an asymptotically 

perfect packing of Ej—d+iiQi "n)- Thus, the DEP, which is the ratio between the total number of decodable 
subspaces and \Er-d+i{(l, 'n)\, becomes zero. The case where r = ^ and v = ^ + d— \ can be explained 
in a similar fashion. 

We now show that liftings of MRD codes, referred to as KK codes, have the highest DEP among all 
liftings up to a scalar in all nontrivial cases. We denote the DEP of a BDD with decoding radius d — 1 
for a KK code in GF(g)'^^("-'^) with minimum rank distance d as Ps.kk(^i v,d — 1). Since the distance 
distribution of a KK code is transparent to the transmitted codeword, we have removed the dependence 
on the transmitted codeword in the DEP for a KK code. 

Corollary 3: Let C be any rank metric code in GF(g)''^ (r < n — r) with minimum rank distance 
d and let C G C. Then if q > 2, r < n — r, or d n — r — 1, Ps{I{C),u, v,d—l) < HqPs_KK{u, v,d—l), 
where H2 = 3.5 and Hq = ^ for q > 2. 

The proof of Corollary [3] is similar to that of Corollary [T] and is hence omitted. Again note that Corollary [3] 
does not hold only for a trivial case. 

If the probability that the received subspace is at subspace distance u from the sent codeword decreases 
rapidly with u, then the overall DEP is dominated hy Pg{C,d+i,v,d—l), where i = I when d—l + v — r 
is even and i = 2 otherwise. Note that u = d+i is the smallest subspace distance that leads to a decoding 
error and ensures u + v — r is even and hence Ns{r,v,u) > 0. Proposition |6] below determines this value, 



and shows that it asymptotically reaches the upper bound in ( 17 1 for KK codes. 

Proposition 6: The DEP of a bounded subspace distance decoder with decoding radius d — 1 for a 
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CDC in Er{q,n) with minimum subspace distance 2d over a symmetric operator channel, provided that 
the received subspace is at subspace distance d + i (i = 1 when d — l + v — ris even and i = 2 otherwise) 
from the sent codeword, is given by 

Ps(C, d + i,v,d-l) = q-id-)(-+^)^^^^Aa{C), (18) 

Id-rl lT+i\ 

where r = '^~'+"~'' . In particular, the DEP for a KK code satisfies Ps,KK(rf+l, v,d)> /fS^- ""T"" (""2'^+- 
when d - 1 + - r is even and Ps.KK(d + 2,v,d) > K2g-'"'t''"'("-2'^+''"'r+'^)-K"-^+i+f ) when 
d— — ris odd. 

The proof of Proposition [6] is similar to that of Proposition [2] and is hence omitted. Proposition [6] 



and ( 17 1 indicate that the upper bounds in ( 17 1 for a KK code are tight, since the lower bounds in 



Proposition |6] and the upper bounds in ( [T7] l differ by a scalar only. 



C. DEP of Bounded Injection Distance Decoders for CDCs 

We now study the DEP of a bounded injection decoder with decoding radius t = \_^^\ for a CDC with 
minimum injection distance d (and hence, minimum subspace distance 2d) over a symmetric operator 
channel. Recall that if the received subspace has dimension v and is at injection distance from the 
transmitted codeword, then it is at subspace distance u = 2fi — \v — r\. Therefore, the subspaces decodable 
by the bounded injection distance decoder are at subspace distance no more than d—l from a codeword 
and can all be decoded by the bounded subspace distance decoder of the same code. In other words, a 
bounded injection distance decoder only decodes a fraction of subspaces decodable by a bounded subspace 
distance decoder. We refine this statement and express the DEP of the bounded injection decoder for a 
code in terms of the DEP for the bounded subspace distance decoder of the same code in Proposition |7] 
below. 

Proposition 7: A BDD with injection decoding radius t corrects the same subspaces as a BDD with 
subspace decoding radius 2t — \v — r\. Therefore, Pi{C, fi, v, t) = Ps{C, 2^ — \v — r\,v, 2t — \v — r\). 

Proof: A subspace of dimension v is at injection distance fi from the sent codeword and at distance 
no more than t than another codeword if and only if it is at subspace distance u from the sent codeword 
and at subspace distance no more than 2t — \v — r\ than another codeword. Therefore, a subspace is 
decodable by a BDD with injection radius t if and only if it is decodable by a BDD with subspace radius 
2t-\v-r\. m 

We hence apply the bounds on the DEP of the bounded subspace distance decoder to derive bounds 
on the DEP of the bounded injection distance decoder of the same code. 
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Proposition 8: Assuming a codeword C of a CDC in Er{q,n) with minimum injection distance d is 
sent over a symmetric operator channel and that the received subspace has dimension v and is at injection 
distance n from C, the DEP of a bounded injection distance decoder with decoding radius t is given by 

^ r 2t-\v~r\ 

n{C,fi,v,t) = — ^VA^(C) V M2fi-\v-r\,s,2w;r,r,v) (19) 

Ns{r,v,2fi- \v - r\) ^ ^ 

for > d — t + \v — r\ and |t; — r| < t, and Pi{C, fi, v,t) = otherwise. Furthermore, if the CDC is the 
lifting of a rank metric code, then 

i L„<?-(*+"?^)("-2'-+*-^)-^(-2*+^) when d = 2t + l, 

\ L,<7-(*+'^)("-'''+*-'^)-'^("-'*+^)-"+^ whend = 2t + 2, 
where L, = K-2^£^<^-!^^ 



For the injection distance, there is no distinction on v like in ( 17) for the subspace metric, because the 
true subspace distance decoding radius is always given by 2t — \v — r\. Proposition [S] indicates that at 
least d — t errors and at least d — t erasures both have to occur for the bounded injection distance decoder 
to decode erroneously. Therefore, the bounded injection distance decoder is more robust to errors. We 
finally remark that analogues of the other results derived for the subspace distance decoder can also 
be derived for the injection distance decoder, however we shall omit them for the sakes of clarity and 
conciseness. 

The DEPs of the bounded subspace and injection decoders for the same code are compared in Figure 
[T] More precisely, we consider a CDC in £'20(91 50) with minimum injection distance d = 9 and error 



correction capability t = 4. the exponents in the bounds in ( [TT] ) and ( 20 1 are shown forr — d+l<v< 



r + d — 1 and r — t<v<r + t, respectively. We clearly see that both decoders have the same radius 
forr — l<t)<r + l, while their performances diverge for larger values of — r|. This illustrates 
the results in Proposition [TJ where it is shown that the bounded injection decoder can be viewed as a 
bounded subspace distance decoder whose decoding radius decreases with \v — r\. These relations are 
further illustrated in Figure |2| where we schematically compare the volumes of decoding spheres around 
codewords for both metrics. 
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12 14 16 18 20 22 24 26 



Fig. 1. DEP of BDDs with tlie subspace distance and the injection distance for n = 50, r = 20, d = 9, t = 4 (log scale). 



Appendix 

A. Proof of Proposition [7] 

In order to prove Proposition [T] we first need a technical lemma. 

Lemma 2: For all U G Eu{q,n), Y.w&E^{q,n) 9r{U, s,W) = ^Jr{u,s,w). 

Proof: For any A € Ea{q,n), the number of matrices in GF(g)'"^" with row space equal to A 
is given by a{m,a). Hence by counting the number of pairs of matrices (R, W), where R, W G 
GF(g)™^", R{R) = U, R{W) = W, and 4(R, W) = s, in two ways, we have a{m,u)g^{W, s,U) = 
a{m, w)g^iU, s, W). Hence EweEUg,n) 9^^^^ s, W) = EweEUq,n) 5r(^, s, U) = Jr('w^, s, u). 

Using (2), we obtain T.W(iE^{q,n)9^i^,s,W) = |^ Jr(u, s, -«;). ■ 

We now give the proof of Proposition [T] whose outline is as follows. We first prove Q by counting the 
number of decodable matrices. Then, deriving an upper bound on the distribution of rank metric codes 
leads to Finally, we use geometric properties of the rank metric and Lemma [2] to obtain the general 
upper bound in (|7]). 

Proof: Let C be a rank metric code in GF(g)'"^'^ with minimum rank distance d. We have Pr(C, U) = 
Yli=o '^(^' -^i where 6(U, s, C) is the number of matrices X at distance s from the code, and 
such that R{X -C) = U. Hence 6{U, s, C) = ^W:dini{w)>d'^wiC)gR{U, s, W) and we obtain ([5|. 

We now give an upper bound on Ayy{C)- We can express any matrix E G GF((7)™^" with row space 
W e E^{q,n) as E = FW, where F G GF(g)™^"' has rank and W G GF(g)"'^" is a fixed matrix 
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dimension 



dimension 



v+2a 



v+a 



v-a 



v-2a 




v+2a 



v+a 



v-a 



v-2a 




(a) in projective space with subspace metric 



(b) in projective space with injection metric 



Fig. 2. Comparison of the ball with subspace radius 2a (solid line) and the ball with injection radius a (dashed line), both 
centered at a subspace with dimension v (we assume Q < v — 2a < v + 2a < n) 



with row space W . Let Cw = {F : FW + C G C}, then \Cw\ = ^iy(C) and Cw forms a constant-rank 
code in GF(g)™^"' with rank w and minimum distance d. Therefore, Ay/{C) < Af^{q,m,w,d,w) < 
a{m, w — d + 1) by (j4]l. The DEP hence satisfies 



t n 



P^iC,U,t) < V Va(m,u;-d+l) V gK{U,s,W) 

VKeE™(<j,n) 



=0 w=d 



^ n - - t 

— — > a(jn,w - d^\)y JAu,s,w), 



w=d 



(21) 



s=0 



where (2]_i follows Lemma |2] 

When C is an MRD code with minimum rank distance d, it can be shown that A\y{C) = M{q, m, w, d, w) 
[^] ^M{q, m, n, d, w) |7, Lemma 3], which leads to (|8 1. 

By definition of a{m,u), 



a{m, w — d + 1) 



a{m,w)^ ^ , , , < K-^q-'^^^'-'^N^iw), (22) 

a[m — w + d — l,d — 1) ^ 
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and hence 

t n 



p,{c,u,t) < J] J]iVRHJR(u,s,^) 

s=Ow=d 

t n 



=0 w=d 



< K-\-^^^-^W^{t), (24) 



where (23l and (24i follow S and 0, respectively. Using K(t) < |12|, we obtain Fl 



B. Proof of Proposition |4] 

Proof: The outline of the proof is as follows. For A G Ea{q,n) and B € Eh{q,n) with (is(^, -B) = 
u + s, we want to count the number of C G Ec{q,n) satisfying ds{A,C) = u, and ds{B,C) = s. We 
first prove that C DAdB = ACiB and C = C n A + C n i?. We can then count the number of possible 
choices for C n A and C H B, which yields the number of possible choices for C. 

We have dim{C n A) = ^^±f^ = a', dim(C n -B) = ^±f^ = b', and dim(A n B) = "+^-"-" = 
a' + b' — c. This implies that Js(m, s, n + s; a, 6, c) = when u + s>a + b, u>a + c, or s>b + c. 
Since 

dim(Cn^nS) = dim(C n ^) + dim(C n B) - dim(C n A + C7 n 
> a' + b' -c = dim{AnB), 

we have CnAn5 = AnBand dim(C n A + C n B) = c, hence C = CnA + CnS. 

Since C n A is a subspace of dimension a' of A which contains An B, there are N^{q, a, a' + b' — 
c, a', c — b') choices for CCiA. Note that N^{q, a,a' + b' — c,a',c — b') = when s > 2a + b — c. Similarly, 
there are Ns{q, b, a' + b' - c, 6', c - a') choices for CnB. Note that Ns{q, b, a' + b' - c, 6', c - a') = 
when u > 2b + a — c. Thus, there are A''s(g, a, a' + b' — c, a', c — b')N^{q, b, a' + b' — c, 6', c — a') = 
^a-a'^-6'+cj [fc-'^'^-fe'+^j choices for C. m 

C. Proof of Proposition [5] 

We first prove a technical lemma. 

Lemma 3: For all r, s, and t < min |r + s, |_|J }, X]d=o ^sif, s, d) < Lqq-f^'^'^'*\ where 4/(r, s, t) = 
t{2n - t) - {r - s){2n - r - 3s) and Lg = K'^ . 
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Proof: By definition of Ns{r, s, d) and we have N^{r, s, d) < K~'^qf('^''''^\ and hence ^^^^ N^ir, s, d) < 

We now prove Proposition |5j whose outline is similar to that of Proposition [T] After proving the 
equality ( [TS] ), we give an upper bound on the distance distribution of liftings of rank metric codes to 
obtain ([16]). We finally use geometric properties of the projective space and Lemma |3] to obtain the 



general upper bound in ( 17 1 



Proof: Given u and v such that Ns{r,v,u) > 0, we have Ps{C,u,v, d — 1) 



D(C,M,i>) 
Ns{r,v,u) 



, where 



D{C, u, v) is the number of decodable subspaces with dimension v and at subspace distance u from 
C. A subspace is decodable if it is within the decoding radius of any codeword. For a codeword C at 
subspace distance 2w from C, there are exactly X]s=o Js{u, s,2w,r,r,v) subspaces with dimension v, 
at distance u from C, and at distance < d — 1 from C by Proposition [5] Summing for all C, we obtain 

Let /(C) be the lifting of a rank metric code C in GF(g)''^ For C G C, {/(D - C) : D e 
C, dR(D, C) = w} is the lifting of a constant-rank code in GF(g)''^("~^') with minimum rank distance at 
least d and rank w, and hence A«,(/(C)) < Af.{q, n — r,r, d, w) < [^] a{n — r,w — d + 1) by (B, which 



leads to (16 1. 



We now prove ( 17 1 when d— 1 + f — ris even, and the other case is similar and its proof is omitted 
We have a{n - r,w - d + 1) < q-(n~r){d~i)qn,^ ["-rj ^^^^^ ^w) = q""" ^ [""^j and hence 

A^iliC)) < (7-("-^^-)(^-i)iVs(r,r, 2w). 



We obtain 



d-l 



P,{I{C),u,v,d-l) = ^A^{I{C))Y^ 



Js{2w,s,u;r,v,r) 
Ns{r, r, 2w) 



< 



< 



w=d s=0 

r d-l 

^-(n-r)(d-i)^^j^(2u;,s,^z;r,^,. 

w=d s=0 

d-l 

s=0 



(25) 

(26) 
(27) 
(28) 



where (26 1, (27 1, and (28 1 follow (12i, (25 1, (Hsl, respectively. We obtain (HTb since I]f=d ^s(f , r, s) < 



L^q\{d-i){2n-d+i)-\(v^r)(2n-v-^r) Lemma g 
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